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1. Introduction 



The random matrix models liave various applications in rather different domains, 
and sometimes language barriers prevents the flow of ideas and knowledge from one field 
to another. For example, such powerful techniques as the conformal field theory (CFT) 
description of the random matrix models and their relation with the integrable hierarchies, 
which were developped extensively by the string theorists in the early 90's, are practically 
unknown to the mesoscopic physicists. This lecture is an attempt to explain the uses of the 
the CFT description in a language, which is acceptable by both communities. We therefore 
avoided any "physical" interpretation and concentrated on the method as such. The only 
thing the reader is supposed to know are the basics of the two-dimensional conformal field 
theory. 

o o o 

The statistical ensembles of random matrices of large size (matrix models) have been 
introduced in 1951 by Wigner in order to analyze the spectral properties of complicated 
systems with chaotic behavior In this approach the Hamiltonian of a chaotic system 
is considered as a large matrix with random entries. Consequently, the analytical studies 
of random matrix ensembles carried out in the next 25 years (see the Mehta's book |]2|) 
were oriented to the calculation of the spectral correlation functions or joint eigenvalue 
probabilities. 

If M is an X hermitian random matrix, then the spectral correlation function 
p(Ai, An) is defined as the probability density that Ai, are eigenvalues of M 

p(Ai, A.) = ^^^^^ (jl 5{M - A^y (1.1) 

All spectral correlation functions are expressed as determinants of a single kernel K{\, /j.) 
(the spectral kernel). The spectral kernel can be evaluated by the method of orthogonal 
polynomials . Its large asymptotics is characterized by the interposition of a smooth 
behavior and fast oscillations with wavelength ~ 1/A^ (in a scale where the total range of 
the spectrum is kept finite). The smooth large distance behavior depends on the concrete 
form of the matrix potential. On the contrary, the microscopic behavior characterized by 
oscillations depends only on the symmetry group (the unitary group in the case) and fall 
into several universality classes. It is the microscopic behavior of the spectral correlations, 
which is interesting from the point of view of applications to chaotic systems. A review of 
the latest developments in this direction, in particular these related to the recent monte 
Carlo data in lattice QCD, can be found in [Q. 

o o o 

The discovery by 't Hooft of the 1/A^ expansion Q gave a meaning of the smooth part 
of the spectral correlators and opened the possibility of using random matrix models to 
solve various combinatorial problems, the simplest of which is the enumeration of planar 
graphs PJ§]. Here we can mention the exact solution of various statistical models defined 
on random surfaces J^jH^il? the matrix formulation of the 2d quantum gravity ||TO|JTlJ| and 



some more difficult combinatorial problems as the enumeration of the "dually weighted" 
planar graphs |T^ and the branched coverings of a compact Riemann surface |]1^ . 
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In this kind of problems the solution is encoded in the expansion of the loop 
correlation functions, which are the correlation functions of the collective field variable 



^(^)=*M— ^j- (1-2) 



In the large limit the correlation functions of the resolvent (|1.2|) are meromorphic func- 
tions with cuts along the intervals where the spectral density is nonzero. The discontinuity 
along the cuts gives the smooth part of the joint eigenvalue probabilities Indeed, we 

have the evident relation 

{W{z,) . . . W{z^))^ = I f;---f" , . p(Ai, AO (1.3) 

where ( )c means connected correlator. 

The first exact results in the large limit were obtained by direct application of 
the saddle point method P, but later it was recognized that a more powerful method 



is provided by the so called loop equations whose iterative solution allows one to 
reconstruct order by order the expansion. The most efficient iterative procedure 
proved to be the "moment's description" [^, which allowed to calculate the free energy 
and loop correlators for an arbitrary potential up to terms. 

o o o 



In the early 90's, after the publication of the seminal papers [|TT|], the two resolu 



tion techniques in random matrix models (orthogonal polynomials and loop equations) 
developed rapidly and were recognized as particular cases of well developed mathematical 
methods. 

The method of orthogonal polynomials was reformulated in terms of the theory of 
r functions of integrable hierarchies (see, for example, the review article [1^). On the 



other hand it was observed that the loop equations generate a representation of (half of) 
the Virasoro algebra and are therefore equivalent to the requirement that the theory is 



conformal invariant [|T^,|T^. This allowed to describe a class of large N matrix models 
near criticality in terms of the Hilbert space of twisted bosonic fields, and apply the well 
developed formalism of the conformal field theory. Finally, it was observed that the two 
approaches are closely related since the r functions associated with the matrix models 
can be formulated as fermionic theories with conformal invariance (The approach 



based on integrable hierarchies is however less general since it works only for the unitary 
ensemble of random matrices.) 

o o o 

It happened that the CFT techniques were developed exclusively from the point of 
view of their application to noncritical string theories, and therefore are known only to 
a relatively narrow circle of physicists. The aim of this lecture is to demonstrate the 
technical and conceptual advantages of the conformal field theory description in a more 
general setting. In particular, we will reproduce, using the CFT formalism, the microscopic 
oscillations of the spectral kernel and correlators. The CFT formalism allows to understand 
better the origin of the observed universality at mesoscopic and microscopic scales. 
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We will restrict ourselves to the simplest example of the hermitian one matrix model 
with arbitrary potential and will only briefly consider the generalizations at the end. First 
we will establish the equivalence of the unitary matrix ensemble with a system of two- 
dimensional free chiral fermions. Then we will construct the collective fleld theory, which 
is obtained from the fermionic system by the two-dimensional bosonization rules. The fleld 
posesses a nontrivial expectation value associated with a hyperelliptic Riemann surface. 
The latter is described, from the point of view of CFT, by a collection of twist operators 
associated with its branch points, with interaction concentrated at the branch points. We 
will demonstrate the power of the CFT description by performing some actual calculations. 
We will obtain the quasiclassical expressions for the correlator of the resolvent W{z), the 
joint eigenvalue probabilities and the spectral kernel. Finally we will briefly mention how 
the formalism can be generalized to the case of a chain of coupled random matrices in 
presence of external matrix flelds. 



2. The hermitian one matrix model. Loop equations as Virasoro constraints 

The partition function of the unitary ensembleis deflned as the integral 

Zn[V]= J c/Me-T--^(^) (2.1) 



where dM denotes the translational invariant measure in the space N x N hermitian 

V 

N 



matrices M = {Mij}f-^. We normalize the measure as 



dM 



1 rfM^ 2 dRe(M,,) rfIm(M,,) (2.2) 



Vol[t/(iV)] 11 2 



k<j 



where Yol[U{N)] = Y[k=i k\ volume of the unitary group. We assume that the 

potential is an entire function 

oo 

V{M) = -J2tnM^ (2.3) 

n=0 

which is the natural choice for the ensemble of hermitian matrices. 

The integrand depends on the matrix variable M only through its eigenvalues 
Ai, Atv and the integral ( |2.1| ) is actually reduced to 

1 r ^ 

2n[V] = — jUdX, e-^(^») HiX, - X,f. (2.4) 

1=1 i<j 

The loop equations represent an inflnite set of identities satisfled by the correlation 
functions of the collective fleld variable 

N 

i=l * 
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The derivation goes as follows. From the translational invariance of the integration measure 
dXi we find 

J2i^+^J2v^ + J2^'->^r')^) =0 (2.6) 

^^^^^ j^^)^'-^^ ^0 ^^-A./iv,w 

where ( )N,{t} means the average with respect to this partition function. Using the 
identity 

V - + 2 V — - = V — — (2 7) 

we write (^]^) as 



N 



/-^ / N,{t} 

i=l * n>0 ^ ' 

Assume that the point z is very far from the origin. Then the sum in the second term 
can be viewed as the result of a contour integration along a big contour C circling all the 
eigenvalues Ai, Aat but not the point z: 

-^■-^(t{z')) =0 (2.9) 
c2mz-z'\ ^ '/N,{t} ^ ' 

where 

T{z) = \[d^{z)]\ Hz) = ^Y.^^z--V2tT\og{j^). (2.10) 



n>0 



The insertion of the operator trM"^ = can be realized by taking a partial derivative 

with respect to the coupling tn- This allows to represent the operator (|2.10| ) as 

T{z) = -[dHz)f. Hz) = ^Y.^nZ- + V2N\nz + V2Y,'-—^ (2.11) 

^ n>0 n>0 

and write the loop equations as linear differential constraints 

L„-Zjv[t]=0 (n>-l) (2.12) 

where 

(We used the notation Z]s[[t] instead of Zn[V] since a concrete representation of the the 
potential V{z) by the set {in}$^Lo °f coordinates (coupling constants) is chosen.) The 
operators L„ satisfies the Virasoro algebra 

[L^, L„] = (m - n)Lm+n (2.14) 
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which imphes that the this matrix integral reahzes a representation of the conformal group. 
In writing the Virasoro constraints we introduced the derivative with respect to the trivial 
variable to, which satisfies 

—ZN[t] = NZN[t]. (2.15) 

OlQ 

It turns to be much simpler to use apparatus of the conformal field theory instead 



of solving directly the loop equations (2.12). The formal solution of these equations is 



given in terms of a chiral Dirac fermion or, through bosonization, of a chiral bosonic 
field with Liouville-like interaction, which is by its definition conformal invariant. Such a 



construction has been proposed in |1S,19 



In the next sections we consider the CFT interpretation of the matrix model from 



more pragmatic point of view than in refs. [18. T^]. Our aim is explore the large limit of 



the operator solution of the loop equations through the conformal field theory and repro- 
duce the 1/A^ expansion for the free energy and the correlation functions. The bosonic field 
can be considered as the collective field describing the Dyson gas of eigenvalues. The 1/A^ 
expansion can be constructed in terms of a free bosonic field defined on a Riemann sur- 
face determined by the classical background (the spectral density in the ther mo dynamical 
limit). The spectral correlations can be determined directly from the correlation functions 
of this bosonic field. 



3. Fermionic representation 

3.1. Fock space representation of the hermitian matrix integral in terms of Dirac fermions 



The model ( p.l| ) was first solved by the so called method of orthogonal polynomials 
2^,^. The method is based on the possibility to consider the matrix model as a system 



of free fermions. Later the fermionic partition function was identified as a r-function of 
the KP hierarchy [^. (For a review of the theory of tau- functions see ||22|| .) 



There are different fermionic representations of the matrix integral depending on the 
choice of the fermionic wave functions. Below we will describe the one that is most natural 
from the point of view of the conformal symmetry. Let us introduce the chiral fermions 

^^^) = (£mI' ^hz) = {r^'Kz):r^^Kz)]. (3.1) 



Choosing a local coordinate 1/^ at 2; = oo we expand 

,l:M{z)= ij^r^^z-"--^, V*'"'(^)= Yl '^r^^^z-''-'^. (3.2) 

where the fermion modes satisfy the anticommutation relations 

[lpi''\lpl^''^]j^ = 5a,h5r+r',0- (3.3) 
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The left and right vacuum states of given charge I = {l' \f ^) are defined by 

^ 1 ^ (3.4) 

=^;i^^\l) = (r >/(-)) 

and the corresponding normal ordering is denoted by : :. The action of the chiral fermions 
IS and the energy-momentum tensor is T = hldip^ip — ip^ dip]. 

The fermion bilinears 



(3.5) 



jia) (^^) ^*(a) {^z) : (a = 1, 2) 

J+{z) = tp*(^^tp(^^{z), J_{z) = tp*(^^ij^^^{z) 

satisfy the ^(2) current algebra. 

The currents J''°-'>{z) = '^^Jn^z~'^~^, Jn(a) = ^ i ^ '>pr''°'^ '■ commute with 

the fermionic fields as 



/ (a) 

[•Jn 5 •Jm J — I ''^n+m,0- 



[^^"\ J„] = [J<f', ^:^"^] = £Lt (3.6) 



We will also denote the su(2) current components by 

E+{z) = J+{z), 



E_{z) = J.{z), (3.7) 



and the w(l) current by 



Hiz) =l[ja)(^)- J(^)(^)] 



#(z) = ^[J(i)(^) + J(^)(^)]. (3.8) 



The SU2 currents are expressed in terms of the Pauli matrices as 

E± =: i/j*a±il' : . 



(3.9) 



Now we are ready to write the Fock-space representation of the matrix partition 
function ( |2.1| ) in terms of the fermionic fields (|3.1| ) . For this purpose we define the " Hamil- 
tonian" 



H[V] = -j^y{z)H{z) 



(3.10) 

(2) > 



2 

n>0 
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and the "screening operator" 



/oo 
dX E+{X). 
- oo 



(3.11) 



Statement: 

The scalar product 

ZN[V] = {N\e^^^^ e^+\0), (3.12) 
where |0) = |0, 0) and (A^| = (A^, -A^|, is equal to the integral (|]|). 

Proof: 

First we expand in fermion multilinears and move the operator e^'^^ to the right 
using the formulas 

until it hits the right vacuum. The result is0 



(3.13) 



/oo 

rfA7/'*w(A)7/'(2)(A) e-^(^). (3.14) 
-oo 

Then, expanding the exponent in modes 

Q+=ll /dAA'^+'^'-V-';Vi'^e-^(^) (3.15) 

r,r'>0 

and using (|3.3D and (p.4|), we reproduce the original integral. 

The representation (|3.12| ) of the partition function implies that the latter is a r- 
function of the KP hierarchy (this fact has been established in |2^]) and therefore satisfies 
an infinite tower of differential equations, the first of which is the KP equation. 



^ Here we assume that the potential grows sufficiently fast at A ^ ±oo and therefore the 
integral is finite if the fermionic operator is regular at 2; = 0. This is indeed the case since all ipr 
with r > are annihilated by the right vacuum. 
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3.2. Relation to the method of orthogonal polynomials 



The method of orthogonal polynomials consists, in this setting, in choosing a basis 
diagonalizing the quadratic form (|3.15| ) 

J2z^-hy^^ = J2Pn{z)bT\ J^z'^-^i^^lO) = 5^P.(z)6^^' (3.16) 

r>| ri>0 r>i ri>0 

where the polynomials Pn{z) = z"^ + ... satisfy the orthogonal relations 

/oo 
dX P„(A)P^(A)e-^(^) = /in+i5^,n. (3.17) 
-oo 

The operators fen'^' and 6n ' are linear combinations of V'-r' ^-r respectively: 

k ^ k ^ 

where the matrix B is triangular and with with unitary diagonal elements 

S^^ = l, Bnk = Q ioi k>n (3.19) 

and therefore its determinant is equal to one, detS = 1. 
Introducing the set of orthogonal functions 

/„(A)=P^_i(A)e-^^(^), n = l,2,... (3.20) 

such that (/fc, /j) = J dXfk{X)fj{X) = Skjhk, we represent the quadratic form (|3.15|) as 



oo 



m,n=l n=l 



and the partition function equals 

Zn[V] = -^{detBf det,,^(A,/,) = h,h2...hN. (3.22) 

3. 3. Spectral kernel and joint eigenvalue distributions 

A complete set of observables in the matrix model is given by the joint distribution 
probabilities for n eigenvalues (1 < n < A^) 

p(Ai, A.) = ^^^^( n ^i>^k - m). (3.23) 

k=i 



The density probabilities ( |3.23| ) are expressed as expectation values of fermionic bi- 
linears. Indeed, keeping the only relevant term in the expansion e^^ = . . . + /N\ + . . . 
we write the partition function as 



Zn[V] 



dXi...dXr 



N 



k=l 



o)= j dXi...dXn p(Ai,...,An) (3.24) 



from where the representation of the probability measure p(Ai, . . . , Aat) as the expectation 
value of operators E^. Integrating with respect to — n of the A's gives 



p(Ai, A, 



{N-n)\ 

Wi 



(A^-n)! 
iV! 



l[r''K>^kW{>^k)) (3.25) 



where by definition {O) = ^|^'|eg[v^o+"|o°^ • particular, the spectral density is the 
pectation value of the fermionic current 



ex- 



p{x) = j^{r^^^{x)i^^^Kx)). 



(3.26) 



the two-point correlators pO[| 



Since we have a system of free fermions, the average is equal to the determinant of 

p(Ai, A„) = ^ ^• det^xn^(A^, A,) (3.27) 
where 

K{X, A') = (A)V'(^) (A')) . (3.28) 

The two-point function (|3.28|) is called spectral kernel and is expressed in terms of the 
orthogonal functions (|3.2C1| ) as 



N 



K{XA') = Y1 



fniX)UX'] 
^ (/n? /n) 



(3.29) 



Another useful representation of p(Ai, A„) is as the correlators of the Cartan current 
H. One finds, using the operator product expansion for the currents H and E^, 



H{z)e^+\0) 



dX 

z-X 



E+{X)e^' 



This relation generates a set of Ward identities the simplest of which is 

' dX 



mz)) 



z-X 



(E+iX)). 



(3.30) 



(3.31) 



This equation together with ( |3.26| ) gives the the representation of the resolvent VF(z) = 

-^{H{z)) as an integral of the spectral density W{z) = Therefore all the 

information about the spectral correlators is also contained in the discontinuities of the 
correlators of H{z). 
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3.4- Virasoro constraints and Hirota equations 

We give here another derivation of the Virasoro constraints, using the Fock space 
representation of the partition function. The level-1 u{2) current algebra satisfied by the 
operators (|3.5| ) decouples into u{l) and slt(2)i parts. The corresponding Sugawara energy- 
momentum tensor also decouples to two noninteracting pieces 

- ^ : aV'f"') : = T{z)+f{z), (3.32) 



where 



2 

a=l,2 



T{z) = I : {2H'^ + E+E_ + E_E+) : = : H'^{z 



6 ^ ^ ' ' ' (3.33) 

Tiz) = : H\z) : 

and H and H are given by ( P77| ) and ( P78| ) respectively. Both pieces commute with the 
operator e^+, which means that theory is conformal invariant. Only the operator T(z) — 
'Yj^ LnZ~^~'^ produces nontrivial Ward identities. For any n > —1, 

[L^,Q+\= j dX[Lr,,E+iX)] = J dX—(^X^+^E+{X)y (3.34) 

Since our potential diverges as A ^ ±oo, the boundary terms of can be neglected, and we 
have 

(A^|e^[^] T{z) e^+|0) = {regular function at z = 0}. (3.35) 

On the other hand, by commuting L„, n > — 1 with e^^^^ we find the nontrivial linear 
differential equations ( |2.12| ) known as Virasoro constraints. 

The partition function is by construction a r-function of the KP hierarchy [|18[ , Zy [t] = 
T~N[t], where the "times" tn are the coefficients in the expansion of the potential V{z) = 
— Yl'^=i inz"^- It is a particular case of the "general solution" of the KP hierarchy presented 
in ||2^ and satisfies the bilinear Hirota equations 

rf^ JJe^*""*")^" (det(z-M)) (det — ^) =0 (A^' < A^). 

n \ / N,{t} \ z-M / 

or 

I dz {y+{z) ■ ZN[t] \ (v_(z) ■ ZN'[t']) =0 [N' < N) (3.36) 

where the vertex operators V± {z) = exp ± Y.n=o tnZ"^^ exp =f In i =f ^^=1 ^ sfc) 
are understood as formal series in tn and They produce insertions of fermionic 

operators in the scalar product. Namely, V+{z)ZN[t] = {N\'i{}^^^^,^e-^^^^*w{z)G\0), 

V-{z)ZN[t] = {N\'i/j2]lie~^^^^(^^ {z)G\0) . The direct derivation of the Hirota equations 
from the matrix integral is given in ||24|| . 
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The "susceptibility" u[t] = 2^\ogZN satisfies an infinite hierarchy of differential 

equations with respect to the "times" tn the first of which is the Kadomtsev-Petviashvili 
(KP) equation 

d'^u d [ du du d^u~ 
4- h6w- h 



dtl ' dti 



0. (3.37) 



dh dti dtp 

One can check directly that it is satisfied by the asymptotic expansion around the gaussian 
point 

The Hirota equations themselves are not sufficient to determine uniquely the partition 
functions -Z^Ar[t]. The missing information is supplied by the first two Virasoro constraints 
Lq^-iZ = describing the reaction of the integral ( |2.1| ) to translations and rescaling of 
the variables, 

cx) oo 

Lo = iV2 + ^nt„a,, L_i=tiN + J2in + l)tn+idn. (3.39) 

n=l n=l 

Let us consider the example of the cubic potential V{z) = —tiz — t2Z^ —Hz^ The 
partition function depends, due to the scaling laws (|3.39|) only on a single combination of 
the three couplings: 

_ tit2 I +2 _ oj. ^ 

and the susceptibility depends on the three couplings as 

w(ti,t2,t3) = 2dl\ogZN = tf^^fix). (3.40) 
The KP equation reduces to an ordinary differential equation 

a / „ „ „, „ „ „, 3 



^(/ + 2/'x-9//'-|r(x))=0. 



4. Bosonic representation 

In most of the applications associated with the large limit we are interested not in 
the exact expressions in terms of determinants of strongly oscillating fermionic eigenfunc- 
tions (typically they have nodes), but in the smooth effective expressions. In this sense 
it is more advantageous to develop the theory of the collective excitations of fermions for 
which the small parameter plays the role of a Planck constant and whose ground 

state is the equilibrium state of the Dyson gas. The free boson representation is very useful 
to study the large limit of the theory. 
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By the 2D bosonization formulas we represent the chiral fermions as vertex operators 



where 



:^*(a)^(a) .-Q^(a)^ (a=l,2) (4.1) 



(^(a)(^) =^(a) +p(a)ln^-^^Z-" (4.2) 

are hofomorphic two-component scalar fields with 

[J^\j}:']=n5a,bSn+m,0, \p'''\ Q^'^ = 5a,b (4.3) 

do that their operator product expansion is then ip^"^) {z)ip^'''> {z') ~ da,b^^{^ ~ ^')- The 
normal ordering -.Jn^Jm- is defined as Jn^ Jm for n < m, Jm Jn^ for n > m and 
: qp : = : pq := qp. The boson Fock space is generated by the oscillators with negative 
frequencies applied to the vacuum vector |0) such that 

p«|0) = l, j;j|0) = (n>0). (4.4) 

The left vacuum (0| is similarly defined, with the normalization (0|0) = 0, and the state 
{N\ is constructed as 

(iV| = (0|e^^^''e-^^'''l (4.5) 
It will be convenient to introduce the fields 

associated with the su{2) and the u{l) algebras. The corresponding currents H = 
H = E± = J± are given by 

f (4.7) 

The Hamiltonian and the screening operator are represented as 

H[V] = Y.inHn=^<f ^V{z)d4>{z), (4.8) 

/oo 
dX : e^^(^) : (4.9) 
-oo 
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and the Virasoro generators are realized as follows 

T{z) = LnZ-"--' = ^ : dcP{z)d<P{z) : (4.10) 

n 

The C/(A^)-invariant correlation functions of in the matrix model are obtained through 
the identification 

Tr-^j^ = V2 d(j)^{z), det(z-M) =: e^'^+^") : (4.11) 

where (j)+{z) is the singular at 2; = part of the operator (j){z). 

The Virasoro constraints ( p.l2|) follow from the representation of the bosonic field 
(f){z) by the operator 

^ -L ^i„." + + viA,,,^ £_:|^ (4.12) 

V n>0 ^ n>0 

acting on the partition function. The bosonic representation is easily generalized for the 
orthogonal and symplectic matrix integrals, see for example ||26|| . 



5. Quasi-classics 

5.1. Ground state 



In the large limit the bosonic field (j) = '^^^^^'^^^ develops a large (of order A^) vac- 
uum expectation value (pdz) = {4>{z)) which is a solution of classical Virasoro constraints 



Tc = -d(j)1{z) = {regular function at 2; = 0} 



(5.1) 



satisfying the boundary conditions at infinity 

d(t)c{z) l^^, 



(5.2) 



Therefore the classical current Hr is of the form 



H,{z) = -M(z) y(z), y(z) 



\ 



2p 



(5.3) 



k=i 



where M(z) is an entire function of z. 

The classical current has discontinuities along the intervals [a2k-i, 0'2k]- In the case 
of a single cut, the conditions (|5.1| ) and ( |5.3| ) are sufficient for reconstructing the ground 
state. In the case of several cuts the situation is more complicated. The potential V{z) 
then has p minima along the real axis (we assume that it grows at infinity) so that the 
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spectral density p(A) = ^:X^[Hc{X + iO) — Hc{X + iO)] splits into p discontinuous pieces pk- 
Each partition Ni, Np {Ni + ... + Np = N) giving the distribution of the eigenvalues 
between the p minima defines a metastable state. The tunneling events leading to the 
decay of the the metastable state are of order which means that in the large limit 
the metastable states are actually stable. The true ground state is determined by requiring 
that the tunnelling amplitude between two neighboring vacua vanish, which yields p — 1 
additional conditions fixing the filling numbers A^i, . . . , Np. 

The classical field (pc = "^(^*~^'* ~ V^*'^'*) changes its sign after circling a branch point. 
This means that the the expectation values of the two original bosonic fields ip<-'^'> [z) and 
[z) are given by the two branches of the meromorphic function ( ^.3|) . Therefore one can 
speak of a single field (^{z) = {(/?(^' (z), v?'^) (^z)} defined on the two-fold branched covering 
of the complex plane given by the hyper-elliptic Riemann surface = 11^=1 ~ (^k), the 
two sheets of which are sewed along the p cuts [a2k-i, (i2k]- 



A. One-cut solution 



The classical field depends on N and on the potential V(z) directly and through the 
positions ai, a2 of the two branch points. A simple but important observation (made first 
by V. Kazakov) is that the derivative with respect to of the resolvent of the random 
matrix (the expectation value Hc{z) in our interpretation) depends only on the positions 
of the branch points. Indeed, according to the normalization condition 



Ztvi 



the derivative dHc/dN satisfies 



(5.5) 

and therefore behaves as 1/z at z — > oo. The only analytic function of the form (|5.3|) with 
this asymptotics is 

dN y{z) ^(z-ai)iz-a2y 
If we make a conformal transformation z —>■ u with 

dX 



which maps the twice covered z-plane to the cylinder — ztt < Imw < ztt, the branch points 
are replaces by an orbifold structure. The rotation around a branch point acts as 



1 

1 



or 



'-UJ) = -(/)(w), (/)(-w) = 4){u). 
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According to (|5T^), the classical field is related tob u = u{z) by 



The explicit expression for the classical current is 

1 / dz y{z) V'{z')-V'{z) 



^ ' 2 /^^ 27rz y{z') z - z' 



Expanding at 2; = oo and using the asymptotics Hc{z) ~ N/z + ... one finds the conditions 

dz z^V'{z) 



1 27rz y{z) 



-2iV4,i {k = 0, 1) (5.8) 



which allow to determine the positions of the branch points corresponding to given poten- 
tial. 

B. Multi-cut solutions 

The above analysis can be carried out in the general case of a classical solution with an 
arbitrary number of cuts. Assuming that the branch points are ordered as ai < ... < a2p, 
we have 

a2k ]\f^ 

Pk{X)dX = —. 

a2k-l 

or, in terms of the zero modes of the classical current, 

dz 

H,{z)— = Nk {k = l,...,p) (5.9) 
Ak 27rz 

where Ak is the contour circling the cut [a2k-i, 0'2k]- 

In this case there the p normalization conditions ( |5.9| ) which determine the filling 
numbers A'^i -f ... + A''^ = N associated with the minima of the effective potential. The 
derivatives of the current with respect to the filling numbers 

dH^jz) _ dujjjz) 

form a basis of holomorphic abelian differentials dujk{z) associated with the canonical A- 
cycles encircling the cuts [a2k-iiCi2k] and satisfying 

/ dujj{z) = 6kj. (5.11) 
27r« JAk 
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The period matrix of the hypereUiptic surface is given by the integralsl 

^fcj = TT- / diVj (/c = 1, 1, j = 1, (5.12) 

where the contour Bk encircles the points a2k, ^2/0+1, ••• ci2p-i so that 

AkO Bj = 5k,j, (/c, j = 1, 1). (5.13) 

The explicit form of the abelian differentials is 

dz ^ y{z) /^^ 27rz y{z) 

The functions dco/dz behave as 1/z at infinity and are completely determined by the 
positions of the cuts. The explicit expression for the classical current is 

HJ.) = i X Mil (5.15) 

where the contour C encloses all cuts and the point z. Expanding in I/2; and using the 
asymptotics (^3) we find 

dz z^V'(z\ 

^^4^ = -2iV4,p (fc = 0,...,p) (5.16) 

c 27rz y[z) 

while the normalization conditions (^) give, after integrating along the contour A^, 

-9/ n.(A)^=iV, (fc = l,...,p) (5.17) 

where 

1 /■"^'^ da 

n,(A) = -/ ^y(^). (5.18) 

(It is sufficient to add only p — 1 such conditions in order to determine ai, a2p.) 

The quasiclassical problem with several cuts has been considered in p7| , p8| , |29[| . The 
authors of |j2^ find a unique solution (the true vacuum) because they impose conditions 
equivalent to the requirement that the tunelling amplitudes between the different classical 
vacua are zero. 



^ The period matrix is generally overdetermined. The condition that the the period matrix 
describes a Riemann surface (Shottkey problem) is that the corresponding theta function is a 
tau- function of KP (Krichever's conjecture proved by Shiota and Mulase). 
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5.2. Quasi- classical expansion 

The quasiclassical expansion of the free energy jFjv = In^Ar is actually the 
expansion 

CJO 

J^n[V] = J]A^2-2gjr(9)[y/jv] + nonperturbative terms. (5.19) 

g=0 

It is also called genus expansion because the term jF^f ) represents the sum of all connected 
Feynman diagrams of genus g. 

The leading term {0{N'^)) is simply the action for for the classical bosonic field 



Ar2jr(o) = jd\d\'p{\)p{\') ln(A - A') - j d\p{\)V{\) 

= 1-1 d^zdmc - ^/ ^yi^)duz). 



(5.20) 



The next term (C(l)) is produced by the gaussian fiuctuations around the classical 
solution (f)c ■ It is most easily calculated by considering the Riemann surface as the complex 
plane with two twist operators associated with its branch points. The notion of a twist 
operator has been introduced by Al. Zamolodchikov (see also ^^). The twist oper- 
ators, which are conformal operators with dimension 1/16, can be represented as vertex 
operators 

a+(z) = :e^[^'^'(^)-^'^'(^)]: = : e^^^^^^ : (^ = a^.-i), 
a.{z) = :e^[^'^'W-^<^'(^)]: = : e'^^^^^) : {z = a^k). ^^'^^^ 
The twist operators are not conformal invariant and have to be dressed by the modes 



of the bosonic field in order to obtain translational invariant invariant star operators ||3^ . 
One can construct perturbatively the star operators S{ak) out of the modes of the twisted 
bosonic field (the descendants of the twist operator a±(afc)), by requiring that the singular 
terms with their OPE with the energy-momentum tensor vanish. Then the free energy can 
be evaluated as the logarithm of the correlation function of the 2p star operators. 

Let us find the explicit expression for the 0{1) term JF^^). Up to corrections, 
the star operators read 



1 



S{a2k-i) = [p{a2k-i)] 24cr+(a2fc-i), S{a2k) = [p{a2k)] ^^(y-{a2k) {k = l,...,p) 

where ^(a) ~ <f^dz{z — a)~^/'^d(j){z) is the coefficient in front of the power {z — a)~^/^ in 
the mode expansion of d(f){z) at the point z = a. Using the explicit expression for the 
correlation function of the twist operators we find 

. . r 1 



1 1 / 

jF^i) = log(5'(ai) . ..S{a2p)) = -—^\nfx{ak) + - log f detfcj 

k=i ^ 



a2k-l — CL2j 



This expression generalizes the one-cut solution found in [|T5[ by solving directly the loop 
equations. 

The evaluation in this way of the next term JF^^) requires more work, but is still 
much more simple than the direct iteration of the loop equations. (The result of was 
obtained with the help of Mathematica, and I reproduced all the coefficients, except the 
first one, using only my pen.) 
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6. Two-point correlators for one-cut background 



6.1. The two-point correlator of the resolvent 

The connected correlation function of the resolvent 



M z' -M 
is related to that of the current H = -^d(j)hy 

W{z,z') = l{d<P4z)d<P4z')) 
= I ( {dcP{z)d<P{z')) 



2 V {z- z 



i\2 



(6.1) 



(We remind that ^+ is the singular at ^ = part of the field (f).) Therefore, in order to 
evaluate W{z., z'), we need the expression for the two-point correlator of the gaussian field 
f{z) defined on the Riemann surface with a cut along the interval [ai,a2] = [a,b]. The 
latter can be evaluated as the 4-point function of two current and two twist operators 

{ifizMz')), = {0\^{zMz')a4a)a_{b)\0) 

, , . 6.2 

= In {^{z - a)iz' - 6) - ^{z - b){z' - a)) 

This function has a nontrivial monodromy, and the determination of the square roots 
depends on whether the two fields belong to the same or to different sheets the Riemann 
surface. The correlator can be also determined from the condition that it behaves as 
lii{z — z') at short distances and that it has the same analytic structure as the background 
field fc{z). The sum over the four terms gives of course the two-point correlator on the 
plane: 

5^(^(^)(^)^(^')(^')) = ln[C(z)-C(^')], C(^) = iz-a)iz-b). 

a,a' 

We have to calculate the correlator of the field (j){z) = which is antisymmetric 

with respect to the ^2 monodromy. Taking the alternated sum over the four branches of 
(H), we get 

/N\ 1 y^{z-a){z' -b) - ^{z-b){z' - a) 
^{z - a){z' -b) + y/{z - b){z' - a) 

Differentiating twice, we obtain for the current correlator 



^/{z-a)iz'-b) ^ ^(z'-a)iz-b) 

{dHz)dHz')). = _ V(;'-^ i(fz!!) (6.4) 
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6.2. The quasiclassical expression for the spectral kernel 



The quasiclassical evaluation of the spectral kernel can be done as well, but in this 
case it is not so easy to control the approximation. Again we split the bosonic field into a 
classical and quantum parts and consider the quantum part as a free field on the Riemann 
surface. The kernel ( |3.28| ) is expressed in the free field approximation as a correlation 
function of vertex operators 



K{\ij)= (:e 



(6.5) 



As in the case of the resolvent, we first write the quasiclassical expression for the four-point 
function of two vertex and two twist operators 



g{z,z') = (0| : e'^(^) : : e'^^^') : a , (a)a_ (6)|0) 



1 e" 
2' 



<t>c(z)-<t>c(z') 
^2 



Z — Z 



{z-a){z' -h) 
{z-h){z' -a) 



.1/4 



+ 



{z-h){z' -a) 
{z-a){z' -b) 



l/4^ 



(6.6) 



where we have used that {^{z)) = ±-^(j)ciz), the sign depending on whether the field is 
evaluated on the upper or the lower sheet of the Riemann surface. The overall coefficient 
is determined by the behavior aX z ^ z' where we have to reproduce -^dcpdz) = ±iTTp{z). 
The branch of (pdz) has to be chosen so that g{z) decays exponentially at infinity. 

The kernel K(X, fj.) is defined on the real axis and is obtained as the average over the 
four values of this function on both sides of the cut [a, b] 



1 



^(A + zeO,;U + ze'0). 



(6.7) 



(The average should be taken because of the ambiguity due to presence of the cut; this 
ambiguity of Q appears only in the large limit.) One finds explicitly, for A, /U G [a, b] 



sin ^ 



iV2 



+ 



1 ^{\-a){h-ii)-^ih-\)(p-a) ^(X-a)(b-ii) + y/{b-X)(f,-a) 

2n [{X-a){b- X){^-a){b- ^)]^/^ 



(6.^ 



This expression coinsides with the one obtained in and in ||3j,^ by solving the ap- 
propriate loop equations. 

When A, /U G [a, 6] and A — ~ 1/A^, one obtains the well known short distance 
behavior 



K(A,^) 



sin[7rAr(A-^)p(^]) 
7r(A — n) 



(6.9) 



In the case when one of the arguments is outside the eigenvalue interval, the exponent is 
real and negative and the spectral kernel decays rapidly. 
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7. Generalization to chains of random matrices 



The most natural generalization of the CFT construction is given by the ADE matrix 
models which were introduced as a nonperturbative microscopic realization of the 
rational string theories with C < 1. Each one of these models is associated with a rank 
r classical simply laced Lie algebra (that is, of type Ar, Dr, Eqj^s) or its affine extension, 
and represents a system of r coupled random matrices. 

Here we will discuss only the models of the A-series, for which there exists a simple 
fermionic representation. The model associated with Aj. = su{r + 1) represents a chain 
of r Hermitian matrices Ma of size Na x A''^ (a = 1, r), interacting by means of r — 1 
auxiliary gauge- field-like rectangular complex matrix variables (a = 1, r — 1) of size 
X Na+i- In this way the Ma and Aa are associated respectively with the nodes and the 
links of the Dynkin diagram of Ar- The partition function of the matrix chain is given by 
the following integral 

Z^[V] = jwiidMa e-'^^^^"^^^ jW rfA^rfAje-^'-^^^J^^-^'-^J^^^^+i. (7.1) 

a=l ^ d=l 

After integrating with respect to the A-matrices and the angular variables of the M- 
matrices, the partition function reduces to an integral with respect of the eigenvalues 
Xai, i — 1, Na, of the hermitian matrices Ma 

a=l i i<Cj a=l i,j '-^ 

Important remark: the integral with respect to the A-matrices exists only if all eigenvalues 
of the matrices Ma are positive. This can be achieved by an appropriate choice of the 
potential. We will thus assume that the integration is restricted to the positive real axis 

Aa > oi. 

Let C+ be a contour representing the boundary of the half-plane Kez > 0. Using the 
fact that the eigenvalue integration is restricted to the positive real axis, we can write the 
following loop equations for each a 

Wa{zf+l ^.-^W{z')\y^G'^'W,{-z')-d,V'^{z)\) =0. (7.3) 
Jq^Itxiz — z \/N,t 

The representation of this integral as a Fock space expectation value is a gener- 
alization of the su{2) construction of Sect. 3. Now we consider r -\- 1 fermion fields 
"(/^(""'(z), a = l,...,r + 1, whose modes in the expansion (|3.2|) satisfy the anticommuta- 



tion relations [i/jr^Ki^t ']+ = ^rs^ab, so that the bilinears J<-°-'>{z) =: tl;*'-'''' {z)'il)<-°-'> (z) :, 



o + 

Another way to achieve this is to take Ma = BaBa where Ba are complex Na x Na matrices. 
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Jai^) — i^*'-"^ iz)ijj^ \~^) cind J„ Lz) = tp* {z)'i/j'-''^ {—z) generate an algebra re- 



lated to the u{r + 1) current algebrap. As before, the u{l) current J = Jt'^) 
completely decouples so that the Cartan currents of su{r) are given by the differences 



(a + l) 



Let us express the potentials V^iz) as differences V^iz) = V'-°-^{z) — \—z), 
and define the Hamiltonians H(°-'>[V'^°-^] = § ■^Y^°-'^{z)J^°-'>{z). It is easy to see that the 



partition function (|7.2|) can be written similarly to (|3.12|) as 

0=1 a=l 

where the vacuum state of charge / = (f \ ...J^ ^ ') is defined by ( |3.4| ) with f ^ ^ — /(°-) = 

2Na, and the "screening operators" are given by (5+ = dX [J^°-'>{z) — J ^ (~^)]- 

The bosonic representation is obtained according to ( [4.1|) in terms of the r + 1 bosonic 
fields a = 1, r + 1. They are split into r fields associated with the with the su{r) 
and the w(l) parts 

<P^iz)=^"^^iz)-^'^^'\-z) (a = l,...,r), ^(^) = _L^^(^(.)((-)«-i;^). (7.5) 

r + 1 ^-^ 

a 

It is convenient to define another set of fields (j)a by 

r{z) = 2M^) - ^ C^'M-z) (7.6) 

6 

where = (a" ■ a^) is the adjacency matrix of the Dynkin diagram. The fields cp"' 
and (pa are related as the contravariant and covariant components of a vector field in the 
base of the simple roots 5" of the su{r). Note however the refiection z —>■ — 2; in the second 
term of ( fT^ ) . 

If we define in a similar way the covariant components Va of the potential by 

Viz) = 2Va{z) -Y^G'^Wbi-z) (7.7) 

b 

then 



(Pa{z) = ~Va{z)+Wa{z), Wa{z) = Y,^—. (7.8) 

1=1 

The loop equations ( |7.3| ) read, in terms of the contravariant components cp"" 

dz' \ 
-[dr{z')f)=^ (a = l,...,r). (7.9) 

c+ z- z / 



The difference originates in the minus sign in the definition of the currents corresponding 
to the simple roots of u{r + 1). 
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The Riemann surface defined by the classical solution represents an (r + l)-fold covering 
of the complex plane which has in the simplest case two cuts [a, b] and [—6, —a] on the real 
axis (0 < a < 6). The fields f'-"-^ represent the values on the different sheets of a single 
meromorphic field f>{z). The quasiclassical expressions for the correlation functions and 
the free energy are obtained by introducing 2r twist operators 



a 



+ [{-)'^a] =: e^^'[(-)'«] : , a^lH'b] =: e-^^'^")'^] : (7.10) 



associated with the points ±a, ±b. In the limit a ^ 0,b oo the r twist operators located 
at the points b and —b merge into a twist operator at the origin while the rest form a 
twisted boundary condition at infinity. In this case the bosonic description is given by 
the 2n twisted field considered in refs. |T^ . 



8. Another generalization: The hermitian matrix model as an amplitude be- 
tween two GL{N) representations 

Our CFT representation actually holds for a more general class of matrix models, 
containing external matrix sources. The source can be introduced by replacing the left 
vacuum state with an excited state. Consider the integral 



^.(1) ,H(2) = j rfMe-*=^^(^)x,a) m (8.1) 

where XRi^) ci-nd Xjj{2) are the characters of two polynomial representations of the group 
U {N) . The polynomial representation R is labeled by the Young tableau 

Yr = {h >l2> ... > In} 
and the corresponding character is given by a ratio of two determinants 



XJM) 



det A- 



det A 



(z,i = l,...,iV). 



(8.2) 



Therefore the partition function can be written again as an integral over the eigenvalues 



(1) 

n^A. e-^(^^)det(Ar'^'— )det(A-^+^"' 

i=i 



N-j + l 



(8.3) 



The U{N) characters are related to the expectation values of the fermions as 



n.<,(A.-A,: 
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(8.4) 



where r'j^^ = k — ^ +^jv-A;+i ~ 1' '^^^ Fock representation of the matrix integral 

is an evident generahzation of ( |3.12| ) 

ZRa)Ri2)[V] = (i?(i)i?(^)|e^[^]e^+|0) (8.5) 

where 

N 

(i?wi?(^)| = (0|n^%V't) (8-6) 
generalizes the vacuum state |A^). The analog of the diagonalization formula ( |3.22|) is 



T-riV , 

«{i),H{2)L'^J detPfl(i) detP^(2) ^ ^ 



where the matrix Pr labeled by the Young diagram R = {li > I2 > ■■■ > In}^ or equiva- 
lently, by the sequence of half-integers = /c — i + l^-k+i: is made of the coefficients of 
the orthogonal polynomials P^(A) = Yl=i Pu,k^'^~^- P-]^ = Pn-^j- When = 0, 



the determinant of the triangular matrix Pfi(a) equals 1 and (|8.7| ) reduces to the standard 
formula ( p:22|) . 

The quasiclassical expressions for the spectral correlations depend on the nontrivial 
left vacuum only through the classical field (pc and are therefore given again by ( |6.4|) and 
( |6.8|) . The Virasoro constraints are however more complicated because the left vacuum 
state is not annihilated by the negative components of the current. 



9. Concluding remarks 

The CFT representation of the matrix models originates in the SU{N) invariant ma- 
trix integration measure and therefore can be actually constructed for all matrix models. 
In a system of several or even infinitely many interacting random matrices, one can asso- 
ciate with each matrix variable Ma a current Ja{z) = Tr and a vacuum |0a), which 
is invariant with respect to conformal transformations generated by the Sugawara stress- 
energy tensor Ta{x) ~ Jai^)- 1^ the large limit all other matrices interacting with Ma 
can be reduced to an effective matrix mean field, which will plays the role of the potential 
V in our analysis of the hermitian matrix model. The classical Virasoro constraints ( |5.1J ) 
hold again, but the effective potential can have cuts and poles. The Virasoro generators do 
not have in general a simple representation as linear differential operators as the Virasoro 
constraints ( |2.12| ). The classical background is described by a Riemann surface, which is 
in general not a hyperelliptic one. The Hilbert space associated with the matrix variable 
Ma corresponds to a contour separating the cuts of the resolvent Tr^— from the cuts of 
the effective potential. 

The usefulness of the CFT approach consists in the fact that once the classical solu- 
tion faiz) is found, the loop correlators and the spectral kernel are given by the universal 
expressions ( |6.4|) and ( |6.8|) . In this sense observed universality of the short distance be- 
havior of the loop correlators and the spectral kernel is a rather trivial consequence of the 
conformal invariance. It is also possible to construct a universal loop diagram technique 



for the 1/N corrections, one particular case of which is considered in ref. |]3G 
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